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Abstract 

The effects of uniaxial tensile strain on the ultimate performance of a dual-gated graphene nanoribbon field-effect 
transistor (GNR-FET) are studied using a fully analytical model based on effective mass approximation and 
semiclassical ballistic transport. The model incorporates the effects of edge bond relaxation and third nearest 
neighbor (3NN) interaction. To calculate the performance metrics of GNR-FETs, analytical expressions are used for the 
charge density, quantum capacitance, and drain current as functions of both gate and drain voltages. It is found that 
the current under a fixed bias can change several times with applied uniaxial strain and these changes are strongly 
related to strain-induced changes in both band gap and effective mass of the GNR. Intrinsic switching delay time, 
cutoff frequency, and / 0 n// 0 ff rat i° are a l so calculated for various uniaxial strain values. The results indicate that the 
variation in both cutoff frequency and / on // 0 ff rat i° versus applied tensile strain inversely corresponds to that of the 
band gap and effective mass. Although a significant high frequency and switching performance can be achieved by 
uniaxial strain engineering, tradeoff issues should be carefully considered. 
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Background 

Graphene is a promising material for nanoelectronics 
due to its high carrier mobility at room temperature 
and excellent mechanical properties [1,2]. However, the 
on-current-to-off-current ratio of graphene channel field- 
effect transistors (FETs) is very small due to the lack of a 
band gap. As a result, monolayer graphene is not directly 
suitable for digital circuits but is very promising for ana- 
log, high-frequency applications [3]. A sizeable band gap 
can be created by patterning the graphene sheet into a 
nanoribbon using planar technologies such as electron 
beam lithography and etching [4,5]. The band gap of a 
GNR depends on its width and edge orientation. Zigzag- 
edged nanoribbons have a very small gap due to localized 
edge states. No such localized state appears in an armchair 
graphene nanoribbon (AGNR). Son et al. [6] have shown 
that the band gap of an armchair graphene nanoribbon 
(AGNR) arises from both the quantum confinement and 

Correspondence: gskliros@ieee.org 

Department of Aeronautical Sciences, Division of Electronics, Electric Power 
and Telecommunication Engineering, Hellenic Air-Force Academy, Dekeleia 
Attica GR-1 010, Greece 



the edge effects. In the presence of edge bond relaxation, 
all AGNRs are semiconducting with band gaps well sepa- 
rated into three different families N = 3p, N = 3p+ 1, and 
N = 3p + 2, with p an integer, and in each family, the gap 
decreases inversely to the ribbon width [6] . However, the 
band gap of the family N = 3p + 2 is significantly reduced, 
resulting in a close-to-metallic channel. This classification 
has proved very helpful in the study of AGNRs since inves- 
tigating AGNRs of various widths an equivalent behavior 
of ribbons of the same family is revealed. 

Strain has important effects on the electronic proper- 
ties of materials and has been successfully employed in 
the semiconductor technology to improve the mobility of 
FETs [7]. For GNRs, it has been established that the band 
structure can be drastically modified by strain. As a result, 
it has been proposed that strain can be used to design vari- 
ous elements for all-graphene electronics [8]. The effect of 
strain on the electronic structure and transport properties 
of graphene sheets and its ribbons have been studied both 
theoretically [9-11] and experimentally [12-14]. Uniaxial 
strain can be applied by depositing a ribbon of graphene 
on transparent flexible polyethylene terephthalate (PET) 



Spri 



inger 



© 201 4 Kliros; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution 
License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any 
medium, provided the original work is properly cited. 



Kliros Nanoscale Research Letters 201 4, 9:65 
http://www.nanoscalereslett.eom/content/9/1/65 



Page 2 of 1 1 



and stretching the PET in one direction [12]. Moreover, 
local strain can be induced by placing the graphene sheet 
or ribbon on a substrate fabricated with patterns like 
trenches as it has been explored for achieving quantum 
Hall effect [15]. To date, however, no experimental works 
on applying uniaxial strain to narrow GNRs (of sub- 10 nm 
width) have been reported. 

In comparison to a graphene sheet, whose band gap 
remains unaffected even under large strains of about 20%, 
the band gap of GNRs is very sensitive to strain [16]. 
Since shear strain tends to reduce the band gap of AGNRs, 
most studies are concentrated to uniaxial strain. Uniaxial 
strain reduces the overlapping integral of C-C atoms and 
influences the interaction between electrons and nuclei. 
As a result, the energy band structure, especially the 
lowest conduction subbands and the highest valence sub- 
bands should be changed. Recently, the band structure 
and transport properties of strained GNRs have been the- 
oretically explored using tight binding as well as density 
functional first-principles calculations [16-19]. It is found 
that uniaxial strain has little effect on the band structure of 
zigzag GNRs, while the energy gap of AGNRs is modified 
in a periodic way with a zigzag pattern and causes oscil- 
latory transition between semiconducting and metallic 
states. Moreover, the band gaps of different GNR families 
show an opposite linear dependence on the strain which 
offers a way to distinguish the families. Tensile strain of 
more than 1% or compressive strain higher than 2% may 
be used to differentiate between the N = 3p + 1 and N = 
3p + 2 families as their band gap versus strain relation- 
ship have opposite sign in these regions [18,20]. However, 
shear strain has little influence on the band structure of 
AGNRs. On the other hand, neither uniaxial strain nor 
shear strain can open a band gap in zigzag GNRs due to 
the existence of edge states [16]. 

Although several studies have investigated the band 
structure of strained AGNRs, only a few have been 



focused on the performance of strained GNR-FETs 
[21-24]. These studies are based on first-principles quan- 
tum transport calculations and non-equilibrium Greens 
function techniques. It is shown that the I-V characteris- 
tics of GNR-FETs are strongly modified by uniaxial strain, 
and in some cases, under a 10% strain, the current can 
change as much as 400% to 500%. However, the varia- 
tion in current with strain is sample specific [22]. On the 
other hand, although semi-analytical [25] or fully analyt- 
ical models [26] for the I-V characteristics of unstrained 
GNRs-FETs have been proposed, no analytical model of 
GNRs-FETs under strain has been reported. 

In this work, using a fully analytical model, we investi- 
gate the effects of uniaxial tensile strain on the I-V charac- 
teristics and the performance of double-gate GNR-FETs. 
Compared to top-gated GNR-FET, a dual-gated device 
has the advantage of better gate control and it is more 
favorable structure to overcome short channel effects [27] . 
Since significant performance improvement is expected 
for nanodevices in the quantum capacitance limit QCL 
[28], a double-gate AGNR-FET operating close to QCL 
is considered. High frequency and switching performance 
metrics of the device under study, as transcoductance, cut- 
off frequency, switching delay time, and power-delay time 
product are calculated and discussed. 

Methods 

Device model 

Effective mass and band structure 

The modeled GNR-FET has a double-gate structure with 
gate-insulator Hf02 of thickness t ms = 1 nm and rela- 
tive dielectric constant k = 16, as shown schematically 
in Figure la. The channel is taken to be intrinsic, and its 
length is supposed equal to the gate length Lq. The source 
and drain contacts are heavily doped regions with dop- 
ing concentration value of 5 x 10 -3 dopants per carbon 
atom. Since thin and high-/c gate insulator is employed, 
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Figure 1 Schematics of double-gate GNR-FET and the atomic structure of AGNR. (a) Schematics of double-gate GNR FET where a 


semiconducting AGNR is used as channel material, (b) The atomic structure of AGNR. Hydrogen atoms are attached to the edge carbon atoms to 


terminate the dangling bonds. N is defined by counting the number of C-atoms forming a zigzag chain in the transverse direction. 
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we can expect excellent gate control to prevent source- 
drain direct tunneling. Moreover, the quantum capaci- 
tance limit (QCL), where the small quantum capacitance 
dominates the total gate capacitance, can be reached. The 
channel material is assumed to be a single-layer AGNR of 
the family N = 3p + 1, as it is illustrated in Figure lb. It 
is well known that this family of AGNR is semiconduct- 
ing material with promising characteristics for switching 
applications [26]. The edge boundaries are passivated by 
hydrogen atoms. It has been demonstrated that hydrogen 
passivation promotes the transformation of indirect band 
gaps to direct ones resulting in improved carrier mobility 
[19]. Moreover, the edge of the GNR is assumed to be per- 
fect without edge roughness for assessing optimum device 
performance. In what follows, a power supply voltage of 
y DD = 0.5 V and room temperature T = 300 K are used. 

Before dealing with the device performance under 
strain, we consider the effect of uniaxial strain on both 
band gap and effective mass of the AGNR. It has been 
verified that a 3NN tight binding model incorporating 
the edge bond relaxation can accurately predict the band 
structure of GNRs [29] . The 2NN interaction, which only 
shifts the dispersion relation in the energy axis but does 
not change the band structure, can be ignored. Any strain 
applied into the GNR modifies the C-C bonds accordingly. 
As a result, each hopping parameters in the tight-binding 
Hamiltonian matrix of the unstrained GNR is assumed to 
be scaled in Harrisons form [30] t{ = to(di/do) 2 , where d{ 
and do are the C-C bond lengths with and without strain, 
respectively. Following the analysis of [16], where these 
changes are treated as small perturbations, we can express 
the energy dispersion of an AGNR under uniaxial strain in 
the form 



'C,n ' 



(hv n k x ) 2 



with 
and 



c\(Yi + 73) + 2c 2 /i5 cos(nO) 
+ 2c 2 y 3 [c 3 + (1 - c 3 ) cos(2h<9)] 



(1) 



(2) 



- -syic 2 cos(n0) [c\(yi + 73) 
(hv n f = (3a cc ) 2 x \ + 2c 2Y3 (c 3 + (1 - c 3 ) cos(2w0))] 

- 73 [cm + (ci - 1)73 

+ 2c 2 73 (c 3 + (1 - c 3 ) cos(2w0))] , 

(3) 

where 0 = 7t/(N + 1), =b indicates the conduction band 
and valence band, respectively, N is the total number of 
C-atoms in the zigzag direction of the ribbon, n denotes 
the subband index, and Ec, n is the band edge energy of 
the nth subband. The strain parameters are expressed as 
c\ = l+a,C2 = l + ^,c 3 = (y 3 c 2 + Ayi)/y 3 c 2 (N+l) with 



a = -26 + 36 2 and 0 = -(1 - 3v)6/2 + (1 - 3v)V/4, 
where e and v are the strength of uniaxial strain and the 
Poissson ratio, respectively. Negative e value corresponds 
to the compressive strain and positive e value corre- 
sponds to the tensile strain. The first set of conduction and 
valence bands have band index s = — 1. Due to the sym- 
metric band structure of electrons and holes, one obtains 
for the energy gap Eq )H = 2Ec >n . Also, y\ = —3.2eV and 
73 = — 0.3eV refer to the first- and third-nearest neigh- 
bor hopping parameters and Ayi = —0.2 eV is used for 
the correction to y\ due to edge bond relaxation effect. A 
poissons ratio value of 0.165 is used in this study [31]. The 
electron effective mass of each conduction subband can 
be calculated by using the formula 

^ ( a ~*:\ Kx> ) 1 (4) 



dkl 



and at the bottom of the conduction band is given by 

* E C,n 



(5) 



Fi gure 2 illustrates the dependence of band gap Eq )H 
of the GNRs family N = 3p + 1 on the uniaxial ten- 
sile strain 6. As it is seen, in the range of tensile strain 
0% < 6 < 15%, Eg decreases first and then increases lin- 
early. Therefore, there is a turning point, i.e., as the strain 
increases, there is an abrupt reversal in the sign of dE g /de, 
making the curves to display a V shape. The turning point 
moves toward smaller values of strain as the width of the 
AGNR increases. Moreover, the slope of E g (e) is almost 
identical for various N and the peak value decreases with 
increasing N. The above observations are in agrement 
with the main features revealed by using tight-binding 
or first-principles numerical calculations [17,20]. On the 
other hand, Figure 3 shows the variation of effective mass 
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Figure 2 Band gap variation versus uniaxial tensile strain for 
different (3p + 1)-GNRswith indices p = 3,4,5,6. 
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Figure 3 Effective mass variation versus uniaxial tensile strain 
for different (3p + 1)-GNRs with indices p = 3,4, 5,6. 



at the conduction band minimum with strain e. As it is 
clearly seen, m* has similar strain dependence as Eg and a 
linear relation between m* and E g is expected which could 
be correlated to an inverse relationship between mobil- 
ity and band gap [32]. These effective mass variations is 
attributed to the change in the conduction band minimum 
position under various strain values. 

Device performance 

Assuming a ballistic channel, the carriers with +k and —k 
states are in equilibrium with Fermi energies of the source 
(£fs) and the drain (£fd)> respectively, with £fs = £f 
and £fd = E? — qVj). The carrier density inside the 
channel can be obtained by employing the effective-mass 
approximation and integrating the density of states over 
all possible energies [26] 



hid 



k B T 
2nh 2 



XI y™*n [^-1/2(^,5) + F-l/2(rj n ,D)] , 



n>0 
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where Fj is the Fermi-Dirac integral of order / defined as 



POO 



rdx 



(7) 



+ exp[* - (ri/ksT)] 

and r] n>s = (E ¥S - E c ,n) /k^T, rj n>D = (E ¥D - E c , n ) IHT. 

Considering the electrostatics describing the structure, 
the following relation between the gate voltage and Fermi 
energy £p can be obtained [33] 



Vg(Ey) ~ ^FB 



£f qn w (E ¥ ) 

q Mns 



(8) 



where q is the carrier charge, Q ns is the gate-insulator 
capacitance per unit length of the GNR and Vfb denotes 
the flat-band voltage. The value of Vfb depends on the 



work function difference between the metal-gate elec- 
trode and the GNR, and it can be set simply to zero as it is 
discussed in detail in [34]. The gate-insulator capacitance 
can be calculated by the simple expression [35] 



Qns =N G K€ 0 



w 

Ans 



(9) 



where Nq is the number of gates (Nq = 2 in our DG- 
device), k is the relative dielectric constant of the gate 
insulator, t- ms is the gate-insulator thickness and a is a 
dimensionless fitting parameter due to the electrostatic 
edge effect. In our numerical calculation, a value of a = 1 
is adopted following [35]. The gate insulator capacitance 
increases linearly as the GNR width increases because the 
area of the GNR increases proportionally. 

The bias-dependent gate capacitance per unit length C g 
can be modeled as a series combination of insulator capac- 
itance per unit length Q ns and the quantum capacitance 
per unit length Cq, that is, 



Co 



Qns Cq 

Qns + Cq 



(10) 



The quantum capacitance describes the change in chan- 
nel charge due to a given change in gate voltage and can 
be calculated by Cq = # 2 9hid/9£f where q is the electron 
charge and hid is the one-dimensional electron density 
[33]. Using Equation (6) and writing in terms of Fermi 
integrals of order (—3/2), we obtain [26] 



CQ = m .l/ 2 J2V^[ F -y2(rin,s)+F-3/2(rin,D)] . 

(27th 2 k B T) n>0 

(11) 

Following Landauers formula and Natoris ballistic the- 
ory [34,36], the device current is expressed by a product of 
the carrier flux injected to the channel and the transmis- 
sion coefficient which is assumed to be unity at energies 
allowed for propagation along the channel. Contribution 
from the evanescent modes is neglected. Thus, 



Id = \ Y\ r C/s(£) -ME)) dE 



(12) 



n>0 \_ C,n 

where fs,T>{E) are the Fermi-Dirac probabilities defined as 
1 
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After integrating, Equation (12) yields 
' 1 + exp(^,s) ' 



n>0 



1 + exp(^, D ) 



(13) 



(14) 



For a well-designed DG-FET, we can assume that Q ns » 
Cd and Q ns » Cs which corresponds to perfect gate 
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electrostatic control over the channel [28]. Moreover, car- 
rier scattering by ion-impurities and electron-hole pud- 
dle effect [37] are not considered, assuming that such 
effects can be overcome by processing advancements in 
the future. In what follows, a representative AGNR with 
N = 16 is considered. 

Results and discussion 

In this section, we firstly explore the calculated device 
characteristics. Figures 4 and 5 show the transfer Id — Vgs 
and output Id — Vbs characteristics, respectively, in the 
ballistic regime, for the DG AGNR-FET of Figure 1 with 
N = 16, which belongs in the family N = 3p + 1, for 
several increasing values of uniaxial tensile strain from 
1% to 13%. The feasibility of the adopted range of tensile 
strain values can be verified by referring to a previous first- 
principles study [22,23]. As it is seen from the plots, the 
current first increases for strain values before the turning 
point 6 — 7% in the band gap variation (see Figure 2) and 
then starts to decrease for strain values after the turning 
point. Moreover, the characteristics for 6 = 5% are very 
close to that of e = 9%, and the same can be observed 
when comparing the characteristics of 6 =3% with that 
of 6 = 13%. Note that, in each region of strain values 
(region before the turning point and region after the turn- 
ing point), there is an inverse relationship between the 
current and the band gap values. Similar features in the 
current-voltage characteristics have been observed in the 
numerical modeling of [22,23] under uniaxial strain in 
the range 0 < e < 11%. These features could be explained 
by the inverse relationship between mobility and energy 
gap which results to an increase in carriers velocity before 
the turning point and to the reduction in carriers veloc- 
ity after the turning point [32]. It is also worth noting 
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Figure 4 Transfer characteristics /□ - Vgs for various tensile 
strain values. 




Figure 5 Output characteristics Id - ^ ds for various tensile 
strain values. 



that, at the ballistic transport limit without electrostatic 
short channel effects, the characteristics in Figure 5 are 
independent on the channel length. This result is dif- 
ferent from conventional FETs and can be explained by 
the fact that, under purely ballistic conditions (no opti- 
cal phonon nor acoustic phonon scattering), the scattering 
mechanisms that cause the channel resistance to increase 
proportionally to channel length are neglected here. 

Now, we focus on the effect of uniaxial strain on the gate 
capacitance C g and transconductance g m = 3Id / d Vq of 
the device under study. Uniaxial strain changes the den- 
sity of states and hence changes the quantum capacitance 
Cq of the channel which is directly proportional to the 
density of states. As a result, in the quantum capacitance 
limit, uniaxial strain changes considerably the intrinsic 
gate capacitance Cg. Figures 6 and 7 show C g versus gate 
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Figure 6 Gate capacitance C g versus V GS for various tensile 
strains. 
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bias at drain bias Vbs = 0.5 V and C g in the on-state 
(where Vgs = Vds = Vdd) versus strain e, respec- 
tively. We clearly observe the non-monotonicity of the 
Cg — Vg characteristics arising from the non-monotonic 
behavior of the function F-3/2(#) in Equation (11). A com- 
parison of the curves in Figure 6 reveals that the gate 
bias Vg at which C g peaks depends on the applied uni- 
axial strain. More specifically, the peak values of C g are 
decreased and moved toward lower values of Vq as uni- 
axial strain is increased before the turning point and are 
increased and moved toward higher values of Vq as uni- 
axial strain is increased after the turning point. On the 
other hand, Figures 8 and 9 illustrate the effect of uniax- 
ial strain on the transconductance g m which describes the 
devices switching-on behavior. As it is seen, g m increases 
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Figure 8 Transconductance g m versus Vqs for various tensile 
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after threshold almost linearly with Vqs and does not peak 
at a certain gate voltage but gets saturated. Moreover, as 
uniaxial strain increases, g m drastically increases from its 
value in the unstrained-GNR case, becomes maximum 
around the turning point e — 7% and then decreases at 
a rate lower than that of the initial increase. This behav- 
ior follows the changes in carriers velocity with uniaxial 
strain, as explained earlier. 

Next, we can assess the high-frequency performance 
potential of the device under strain. The cut-off frequency 
fj is defined as the frequency at which the current gain 
becomes unity and indicates the maximum frequency at 
which signals can be propagated in the transistor. Once 
both gate capacitance and transconductance are calcu- 
lated, fj can be computed using the quasi-static approxi- 
mation [38,39]. 

gm 



/tLg = 



2nC 0 



|V d =^dd 



(15) 



It should be noted that a rigorous treatment beyond 
quasi-static approximation requires the inclusion of 
capacitive, resistive, and inductive elements in the calcu- 
lation. In Figure 5, the quantity /t^g> where Lq is the 
channel length, as function of Vg, for increasing values 
of uniaxial tensile stain, is depicted. Assuming a chan- 
nel length of less than Lq = 50 nm, fj exceeds the 
THz barrier throughout the bias window, confirming the 
excellent high-frequency potential of GNRs. Furthermore, 
Figures 10 and 11 show the variation of cutoff frequency 
versus gate voltage and strain e (in the on-state), respec- 
tively. We clearly observe that fy increases rapidly until 
the turning point e — 7% and then decreases with lower 
rate for higher strain values (e > 7%). This is a direct con- 
sequence of both transconductance and gate capacitance 
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Figure 10 Dependence of {fjL G ) on V GS for various uniaxial 
strains. The drain voltage is held constant at 0.5 V. 



variations with strain. Therefore, the high-frequency per- 
formance of AGNR-FETs improves with tensile uniaxial 
strain, before the 'turning point' of band gap variation but 
becomes worse after this point. 

Lastly, we study the effect of strain on the switching 
performance of the DG-GNR FET. Figures 12, 13, and 14 
show the dependence of I on , 7 G ff and I on /Io& ratio on the 
uniaxial tensile strain, respectively. As it is clearly seen, 
the variation of both 7 on and I Q ff is opposite to the varia- 
tion of the band gap with strain whereas the ratio I on /I 0 {{ 
changes with strain following the band gap variation. The 
on-current 7 0 n changes almost linearly with strain whereas 
the 7 0 ff and the ratio 7 0 n/^off changes almost exponentially 
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Figure 1 2 Variation of the on-current / on versus uniaxial strain. 



with strain. Note that the corresponding curves are not 
symmetric around the turning point, e.g., although for 
6 = 12%, the GNR band gap returns to its unstrained 
value; the drain current at this stain value does not com- 
pletely return to that of the unstrained GNR. This can 
be explained by the fact that although the band gap has 
returned its unstrained value, the carrier group velocity 
has been modified because, under tensile strain, some 
C-C bonds of the AGNR have been elongated [9]. 
Fi gure 15 shows the 7 on versus I on /I 0 & plots for vari- 
ous strains which provides a useful guide for selecting 
device characteristics that can yield a desirable I on /I 0 ff 
under strain. As it is seen, increased tensile strain before 
the turning point of band gap variation leads to lower 
^on/^off ratio, whereas, increased tensile strain after the 
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various uniaxial strains. 



'turning point' leads to higher / 0 n/7off ratio albeit at lower 
on-current. 

Intrinsic delay time r s is also an important performance 
metric that characterizes the limitations on switching 
speed and AC operation of a transistor. Once the gate 
capacitance is calculated, r s is given by [28] . 



r s = 



QVdd 



(16) 



where the on-current is the drain current at Vg=Vd = 
Vdd- Apparently, the switching delay time r s has simi- 
lar variation as the gate capacitance has with strain, as 
it is depicted in Figure 16. Moreover, as it is seen from 
Figure 17, the switching delay time abruptly decreases 
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with strain before the 'turning point' of band gap varia- 
tion but increases rapidly after this point. We can say that 
switching performance improves with the tensile strain 
that results in smaller band gap whereas degrades with the 
tensile strain that results in a larger band gap. It is worth 
noting that the switching delay time for the unstrained 
case (e = 0%) is found to be z s ~ 23 fs/nm, that is at least 
three times larger than the corresponding delay time in 
uniaxially strained-GNR case. Figures 18 and 19 show the 
switching delay time r s as a function of on-current I on and 
hn/hii ratio, respectively. For digital applications, high 
lon/hii ratio and low switching time delay are required. 
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Figure 18 Switching delay time t s /Z-g versus on current / on for 
various uniaxial strains. 
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Figure 20 Power-delay time product Pt s /Lq versus uniaxial 
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strains. 



However, when the I on /Io& ratio improves with the applied 
tensile strain, the I on and switching performance degrade 
and vice versa. Another key parameter in the switching 
performance of the device is the power-delay product 
Pr s = (VDD^on)^s that represents the energy consumed 
per switching event of the device. Figures 20 and 21 illus- 
trate the dependence o of power-time delay product Pr s 
on strain and on I on /I 0 {{ ratio, respectively, where similar 
behavior to that of switching delay-time can be observed. 

Conclusions 

We investigated the uniaxial tensile strain effects on the 
ultimate performance of a dual-gated AGNR FET, based 
on a fully analytical model. The model incorporates the 



effects of edge bond relaxation and third nearest neigh- 
bor (3NN) interaction as well as thermal broadening. 
We have focused on the AGNRs family N = 3p + 1 
which is suitable for device applications. The strong mod- 
ulation of I-V characteristics due to the changes in the 
strain is directly related to the electronic structure of the 
GNR channel region, which is modified as a result of 
changes in atomic structure under strain. The on-state 
current, gate capacitance, and intrinsic unity gain fre- 
quency are steadily improved for tensile strain less than 
the 'turning point' value of the band gap V-type variation. 
The observed trends are in consistency with the recently 
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reported results based on tight-binding quantum trans- 
port numerical calculations [21-23]. Switching delay times 
improves with the tensile strain that results in smaller 
band gap whereas degrades with the tensile strain that 
results in a larger band gap. However, when the I 0 n/hii 
ratio improves with the applied tensile strain, the 7 0 n and 
switching performance degrade and vice versa. There- 
fore, although a significant performance can be achieved 
by strain engineering, tradeoff issues should be carefully 
considered. 

It is worthy noting that since purely ballistic trans- 
port and negligible parasitic capacitances are assumed, 
our calculations give an upper limit of the device perfor- 
mance metrics. Moreover, when metal-graphene contacts 
are used, the on-current of the ARGN-FET are degraded 
[40] by lowering the voltage drop on the intrinsic part 
of the device by a factor of ^bal/(^bal + 2i? cont ) where 
7?bal is the intrinsic resistance of the channel and R con t 
is the contact resistances. Furthermore, in the presence 
of metal contacts, the cutoff frequency is degraded since 
the traversal time of carriers is significantly enhanced 
[41]. On the other hand, our approach may underesti- 
mate the actual concentration of carriers in the channel, 
especially for large drain and gate biases, when parabolic 
band misses to match the exact dispersion relation. How- 
ever, we believe that the present fully analytical study 
provides an easy way for technology benchmarking and 
performance projection. Our study can be extended to 
compressive strain allowing negative values of uniaxial 
strain e in our model. However, as it has been demon- 
strated [42], narrow GNRs exhibit a maximum asymmetry 
in tensile versus compressive strain induced mechanical 
instability, that is, the critical compressive strain for buck- 
ing is several orders of magnitude smaller than the critical 
tensile strain for fracture. Such a large asymmetry implies 
that strain engineering of GNR-devices is only viable with 
application of tensile strain but difficult with compressive 
strain. 
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